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Tlj" ■ We derive the simplest form of the k{k^) order graviton self-interaction La- 

grangian density £g(£g) in the weak field approximation. With the diver- 

Qv^ ■ genceless condition, de Bonder gauge and some combinatoric techniques, we 

O ; derive the mathematically independent basis expressions for Cg{C'i) composed 

^' of 8(24) terms. By using the basis expressions C^{C'^) is reduced to 6(18) 



terms. 
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1. Introduction 

Whenever one studies quantum gravity processes in the weak field approximation pro- 
cedure |1[], where the metric perturbation h^y in ^f^j, from the flat space metric 77^,^ is iden- 
tifled with the gravitational fleld, one always encounters the difliculties in dealing with the 
complexity of the algebraic work which stems from the complicated nature of the gravi- 
ton. A large number of indices appear in the interaction Lagrangian since a graviton has 
a pair of indices and it can even experience self-interaction. Furthermore, one should carry 
out numerous permutation and symmetrization operations on each index-pair and over all 
momentum-index triplets of the gravitons in the interaction Lagrangian to get the correct 
vertices. This is due to the symmetry property of the index pair of h^j, and the bosonic prop- 
erty of the graviton, respectively. In 1967, DeWitt derived the graviton self-interaction 
Lagrangian up to 3- and 4- graviton vertices. This Lagrangian has been believed to be com- 
posed of the smallest terms but not been proved. Using this Lagrangian, cross-sections for 
the processes involving 3- or 4- graviton vertices have been evaluated [^,Q without further 
simplifying the Lagrangian except for the marginal correction of the sign of a term [Q in 
the 4- graviton vertex function. As shown in section 2, 3- and 4- graviton vertices can be 
obtained from C{ and C^ . 

In this paper, we present the basis expressions for £^ and £^ and the simplest form of £^ 
and £"1 in terms of the basis expressions. C{ is simplifled from 11 to 6 in the number of terms. 
For the case of £g, 28 terms are reduced into only 18 terms. Though the numbers of terms 
are still rather large, they cannot be simplifled any more because they are expressed by the 
basis expressions. Using our new C]^ and £^, we obtained the same differential cross-sections 
for the graviton-scalar((7s — ^ gs) [Q, graviton- fermion((7/ -^ gf) [^], graviton- massive vector 
hosoia{gW — * gW) ^, graviton-photon((77 — + (77) [Q and graviton-graviton(5^(7 -^ gg) |^ 
elastic scatterings with those of the corresponding references. In this paper, we review 
the formal procedure to obtain the graviton self-interaction Lagrangian in the weak fleld 
approximation in section 2. To deal with the combinatoric properties of the typical terms 
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in the Lagrangian conveniently, we introduce our index convention which is also helpful 
in computer-aided calculations in section 3. Sections 4 and 5 are devoted to explain our 
strategy for getting the simplest expressions. In these sections, the scheme of constructing 
basis expressions for Cg and £g and also the explicit forms of bases as well as the simplest 
forms of £^ and C^ are given in these sections. Conclusion is made in section 6. 

2. Analysis of the problem 

In the classical theory of general relativity , the action for pure gravity is given by 

Ig = 2k-^ I d'x^f^Rix), (1) 

where k = \/W2txGn with the Newtonian constant Gn-, R{x) the scalar curvature, and 
g{x) = detg^^, i.e., the graviton self-interaction Lagrangian density is given by 

Cg = 2K-^y/^R. (2) 

In the weak field approximation, one can write the metric tensor g^,^ as the sum of the flat 
space metric r^^jy and a metric perturbation h^i, as 

gtiu = Vfj-iy + i^h^jy. (3) 

Here h^u can be identified with the gravitational field by giving [mass]"^ dimension to the 
parameter k. The curvature scalar R is defined by 

R = g^''R^.u 
= g'^Rtxu, (4) 

where R^^, is the Ricci tensor and the Riemann-Christoffel curvature tensor R^^^ is defined 



by 



pA _ n -pA _ f) pA I pa pA _ -pa pA (r\ 



and the afiine connection F^^ is defined by 



rL = 7;9^"{d^gau + d^gaf, - d^g^^) (6) 



By using the identity g^ gxu = S^, we get the k expansion of g^'^, 

g^"" = rj^"" -kW"" + K^ /i^" /i^ - K^ W' hih''p + ... (7) 

Expanding g and a/^ up to ic" order , we get 

v^ = 1 + ^/^/i + \t^\h^ - 2/1^^ V) + i^^'l^' - QhW^K^ + 8/i^^/i^ V)> (9) 

Z o 4o 

with the definition h = M. Then we can obtain the leading term of Cg as 
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£3 = 2f,-^[K\~-d,hd^h + -9^/i^9^/i^) + Oi^% (10) 



where we use the harnionic(de Bonder) gauge, 

d,K = Id^h, (11) 

and omit total derivatives. Then k expansion of Lagrangian density is given by 



= CI + kC] + k'CI + ... (12) 

'0 



The leading term C is obviously of the simplest form as 

Cl = -\d,hd^h+\d,h-pd^hi. (13) 

C^ and C^ can be arranged into simpler forms by using some conditions explained below. 
According to Refs. 0,^, they are reduced into 11 and 28 terms, respectively. But we have 
obtained still more simplified forms of C^ and £^ and proved that they cannot be reduced 
any more. 

The framework of our proof is as follows: Let £^{i = 1,2) be the space composed of any 
mathematically possible expressions which has the typical form of the terms inside £* . (in 

3 



this paper any upper indices i,j, k, I are related only to the order of k.) To get the simplest 
form of £*, we should find out the basis 5* for S\ We start our search for B^ from a set 
C* which contains all the possible terms inside any expression in S^. Of course C* spans S^ 
and has all the possible permutaions of indices for the typical terms. Though C* has a large 
number of elements, it is both countable and finite. Without any condition and neglecting 
the multiplication commutativity, the numbers of elements of C*'s are 



(2!)4.4!' 

With the symmetry property of /i^,y and de Bonder gauge condition, we can reject a lot of 
duplicate elements from C*. If C* has n* elements in this step, we can express C* as 

C = |C;^, C2, C3, ..., C„i|. (15) 

We can obtain the useful relations from the fact that the total derivative of the Lagrangian 
can be neglected. Since any A has derivative factors, we can rewrite it by partial integration 
neglecting the total derivatives. In general the result of this partial integration is also an 
element of S^. Therefore, we can express each c* as an linear combination of cl.^s{k = 
1, 2, 3, ...,n*) as 

4 = E^}.4, (16) 

k=l 

where d^j^s are pure numbers. If we solve these n*-simultaneous equations, we can find the 
basis. Not all the equations are independent and some c*'s are rewritten by itself under 
partial integration as 

ci = cj. (17) 

If there are m^ non-trivial equations in the n*-simultaneous equations, there exist n* — m* 
free variables in the simultaneous equations. Therefore we can take any n^ — m* elements 
from C* to make a basis set B^. It is a common orthogonalization procedure. Getting the 



basis expressions, we can obtain the simplest form of the Lagrangians by expressing the 
original £^ and £^ in terms of these basis expressions. Let us find out them step by step in 
the following sections. 

3. Convention Summary 

First of all, we introduce our new convention for expressions appearing in £} and C? . 
Through partial integration and omitting total derivative, every term of £} and C? can be 
transformed into the following typical forms, 

C?^ : (9?/i??9?/i??/i?7/i??, (18) 

where ?'s are the unknown indices. There are 4 and 6 pairs of indices in the typical terms of 
£^ and £^, respectively. Here we neglect the difference between the upper and lower indices 
because every pair of indices is composed of 1 upper index and 1 lower one - contraction is 
always carried out between them. To show the combinatoric characteristics of the indices, 
we use a number in place of a Greek character as an index. Here is an explicit example of 
our convention. 

V ^ (12) 
5aV -^ (3(12)) (19) 

Distinct numbers should be given to the different indices. Innermost parenthesis means h??. 
If there is a double parenthesis, outer one means derivative. Every number in the parenthesis 
means the correspondent index. For brevity, double derivative expression like (1(2(34))) is 
written as (12(34)). This convention is very convenient to read and write in the computer- 
aided calculations. Especially, it shows the combinatoric characteristics of each expressions 
much more clearly. 



4. C^g reduction 

For the case of Cg, the general form is (a(/i/2))(f'('^i"^2))(^i'^2)- After listing all possible 
combinations of indices made of 4 pairs of numbers, we can find the projection among the 
expressions by using partial integration and neglecting total derivatives as 

( a{lil2)){b{mim2)){nin2) 
= - ( Iil2){ah{mim2)){nin2) - (/i/2)(fe("T-im2))(a(nin2)) 
= +{h{lil2)){a{mim2)){nin2) - {Iil2){h{mim2)){a{nin2)) 

+ ( hl2){a{mim2)){h{nin2)). (20) 



Throughout this procedure we extensively use the de Bonder gauge condition to get the 
factor such as (11) or (1(22)) maximally. Rejecting the dependent expressions one by one, 
we get the basis expressions which can only be transformed into itself. They are given by 
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1(22)) 
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1(22)) 
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1(22)) 
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1(23)) 


;i(23))( 


1(23)) 
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1(23)) 
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1(23)) 
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(44) = d>'h\d^h''^h^p 
(14) = d^h\d^h''^h^p 
(34) = d^^h^d^h^^^K^ 
(34) = d^h\d''h/Kp 
(44) = d^h'^'^d^K^h^ 
(14) = d'^h'^^d^K^h^p 

(14) = d'^h'^^d^h^h^p 
(34) = d^h^'^d.h/Kp. 



(21) 



The basis is made up of 8 expressions. This number is smaller than that (11) of the old L^ 
1^. If one expresses £^ in terms of this basis, it is composed of even smaller number of terms 



, i.e., 6 terms as 



'-'q ~ o .^ 1 "*" o •' 3 + ^ /s nfe ~^ fy fs- 



(22) 



5. £^ reduction 

The method of £^ reduction is basically the same as that of C]^. The general form of 
a term in C^ is of the form, (a(/i/2))(^(^i"^2))(^i'T'2)(^i^2)) which has one more pair of 
indices with one more /i??. There are the same number of derivative factors which make 
the transformation rule alike. We can find the projection of a typical term to the others by 
using partial integration and neglecting total derivatives as 

{a{lil2)){b{mim2)){nin2){rir2) 
= -{hk){ab{'mim2))inin2)irir2) 

-(/i/2)(6(mim2))(a(nin2))(rir2) 

-(/i/2)(6(mim2))(nin2)(a(rir2)) 
= +{Khk)){a{'rnim2)){nin2){rir2) 

+ (/i/2)(a(mim2))(6(nin2))(rir2) 

+ (/i/2)(a(mim2))(nin2)(6(rir2)) 

-(/i/2)(b(mim2))(a(nin2))(rir2) 

-(/i/2)(6(mim2))(nin2)(a(rir2)). (23) 

In this case, we can find the projections of a typical term to 5 other terms. Rejecting the 
dependent expressions, we get the basis expressions. 
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(15)(35) = a^/i^a"/i'^^V^- 

(35)(45) = a^/i^a^/i"^/i«'^V 
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J 19 
J 20 

^23 
^24 



(1(22); 


(3(14)) 


(35) ( 


(1(22); 


(3(45)) 


(45) ( 


(1(22); 
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(1(23); 


(2(34)) 
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(1(23); 


(2(45)) 


(14)( 



(45) = d^h^d'^h^^K'^hp^ 
(13) = d^'h-'^d'^h^'^h^^h^^ 
(53) = d^^h^d^'h^'^hp^h^^ 
(45) = a^r"9,/i^„/i'^'^/i;3. 
(55) = d^'h^'^d^K^h^h^ 
(45) = d^^h^'^d^K^h^'^hf,^ 
(55) = d^^h'^'^d^h/h^ph^ 
(45) = d^^h^'^d^h/h^'^hf,^ 
(45) = d^'h''''d^h''''Kahf,. 
(35) = 9^/i'^"9^/i'^'^/i,^/i„. 
(55) = d^h-'-d'^K^h^ph^ 
(45) = d^h-'-d'^K^h^'^h^, 
(45) = d^h-'-d^K'^h^^h^, 
(35) = 9'^/i'^"9'3/i//i^^/i„. 
(45) = d'^h-'^d^hjh^'^h^, 
(55) = d^h''''d,hjh^ph\ 
(35) = d'^h-'^d^hf'^h^ph^,. 



(24) 



The basis is made up of 24 expressions. This number is smaller than that (28) of the old £^. 
If one expresses C^ in terms of this basis, it is composed of even smaller number of terms , 
i.e., 18 terms as 



£ 



L f2 , i f2 , 1 .2 _ 1 .2 , 1 .2 _ 1 .2 ,]_f2 _ 1 /.2 _ 1 .2 

32 8 4 2 2 4 16 8 "^ 2 

+/15 ~ Q'^is "I" 2^'^ "" 4''i8 ~'~ 2 ^^ ^ "^ "*" "" ~'~ 2 



(25) 



This reduction in the number of terms in C^ and £^ can not be done any further. The 
simplest form shows her face after 27 year's latency. Though the idea we employed to get 



the results is somewhat simple, the actual procedure are hard to be carried out properly if 
one sticks to manual work only. Throughout our calculation, we extensively used a symbolic 
manipulation package Mathematica Q which enabled us to execute the powerful pattern 
matching programs. 

6. Conclusion 

The simplified forms of the graviton self-interaction Lagrangian densities are derived to 
the K and k^ order and they are proved to be the simplest. The key point is that we first 
construct the basis expressions for the C^ and C^ and express them by this basis. We pre- 
sented these basis sets explicitly and described not only the practical way to get the basis but 
also the formal proof. We also checked our results by calculating the observables of several 
real processes. To check the 3-graviton vertex, we evaluated the differential cross-sections 
for the processes, graviton-scalar((7s -^ gs) [Q, graviton-fermion(5f/ -^ gf) ^, graviton- 
massive vector boson((7W^ — > gW) 0, graviton-photon(5f7 — > (77) [Q, by using our simplest 
C^. The differential cross-section of the process graviton-graviton((75f -^ gg) [Q is also eval- 
uated in the same way to check the 3-and 4-graviton vertices. We obtained the same results 
as those of the previous graviton self-interaction Lagrangian. With the simplified graviton 
self-interaction Lagrangian density, we can save much time and avoid errors in dealing with 
the numerous indices. And our new method to find the basis expressions for the £g and 
C^ can be applied to any order of k with slight modification. We presented a successful 
application of the powerful pattern matching abilities of the symbolic manipulation package 
Mathematica to dealing with the physical problems involving complicated index relations. 
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